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Abstract. In this paper we made a review of some papers about probabilistic 
regulatory networks (PRN), in particular we introduce our concept of homo- 
morphisms of PRN with an example of projection of a regulatory network to 
a smaller one. We apply the model PRN ( or Probabilistic Boolean Network) 
to the immune system, the PRN works with two functions. The model called 
"The B/T-cells interaction" is Boolean, so we are really working with a Prob- 
abilistic Boolean Network. Using Markov Chains we determine the state of 
equilibrium of the immune response. 



1. Introduction 

The biological process can be modeled using different class of models, but in 
most of the applications differential equation models have been selected because 
the entities can have more than two values. Here we use a probabilistic regulatory 
network model, that it is a discrete model with only a finite number of state and 
activities. In particular, we describe the dynamic of the immune response in hu- 
mans, using the Boolean model called B/T-cells [8], and we added probabilities and 
use the dynamic model of Probabilistic Boolean Network, [lOj [TTJ [T2l Q31 [14] . 

The immune system is separated by functionality in two parts: recognition and 
effector functions. A complex immune system has cells and molecules that give us a 
basic defense against bacteria, viruses, fungi, and other pathogenic agents. Possibly, 
it has tens or hundreds of different types of regulatory and effector molecules. So, an 
important role in the study of this class of system plays the reduction of networks, 
for that reason we introduce the concept of projection of one net to another smaller 
for the future applications. A variety of cell types compose the immune system, the 
most important are the lymphocytes. These cells are created in the bone marrow, 
along with all of the other blood cells, and are transported throughout the body via 
the blood stream. Lymphocytes spend considerable time resident in lymphoid organs, 
such as the bone marrow, the thymus, the spleen, and lymph nodes. Lymphocytes 
are subdivided in two classes: B-cells and T- cells, see Perelson [5]. B lymphocytes 
secrete antibodies, and the main function of T-cells is the interacting with other 
cells. Helper T cells act through the secretion of lymphokines, they made possible 
to transform the B-cells into an antibody-secreting state. Helper T-cells are the 
cells that are predominantly infected by the human immunodeficiency virus, and 
plays a major role in AIDS. Cytotoxic T-cells, are responsible for killing virally 
infected cells and cells that look like abnormal, such as some tumor cells. 

In this paper we study the dynamic of the B/T-cell model giving by Kaufman, 
Urbain, and Thomas, n 1985. In this model they use functions, in which a boolean 
parameter appears, and they obtained the steady state of the system, in our case 
we describe in a more complex way the immune system using probabilities for the 
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two possible functions of activity. Our model can be changed for another more 
complex if we consider a three values model instead of a Boolean model 

2. Preliminaries and projection 

In this section we introduce the mathematical background of the model Proba- 
bilistic Regulatory Network, and the concept of projection using and example. For 
the complete mathematical background we suggest to see, [3113 ?]. Here, we give a 
method that permit us to build regulatory networks with probabilities assigned to 
its functions. We use an algorithm for understanding the concept of Probabilistic 
Regulatory Network. 

2.1. Algorithm. Input: 

1. n — number of entities in the network under studying, for example 100 genes, 
and the set of values for each entity, that we denote by k a . 

2. A set of relations {m a ,b} taking 1 if the entity a is related to the entity 6, and 
otherwise. 

3. A set of finite families of states in the network which gives the time series data 
for one, two or more update functions, A\ — {(a^ 1; . . . , cq n _ x , a} n )\l < i < m x }, 

and A s = {(a| il ,...,af jn _ 1 ,o? n )|l < i < m s }. 

4. A set of values C = {c\, . . . , c s } with s probabilities obtained in some way by 
the experiment or by the time series data. That is c\ + • • • + c s = 1, and Cj 6 [0, 1]. 

(Alml) Creation the low level graph T: 

1. n = 3 , and m\ 2 = 1 , m 23 = 1 , y TO13 = 0, then our net is very simple 
and it is the following: 

• 3 

r I , 

1 • .2 

(Alm2) We define the set where the functions are acting, in our case, they are 
boolean, that is k\ — {0, 1} = k 2 = &3 = ^2, considering k = k± x /c 2 x &3 = 
{0, l} 3 . The two sequential states are the following: 

Time series data 1 A 1 = {(0, 1, 0), (1, 1, 1), (1, 1, 0), (1, 1, 1)}. 
Time series data 2i 2 = {(0, 1, 0), (0, 1, 1), (0, 1, 0), (0, 1, 1)}. 

We obtain two different update functions: 

fi(xi,x 2 ,x 3 ) = {l 1 x 2 ,x 2 (x 3 + 1)), and / 2 (xi, x 2 , x 3 ) = (xi,x 2 ,x 2 (x 3 + 1)). 

(Alm3) We assign the following probabilities to each update function: 2/3 to the 

function /1, and 1/3 for the function f 2 . 
(Alm4) We construct the high level digraph, that is the following in this example 

Output: 

States Space 

o(Tao) (i7?T) (oTiTo) '3 (oTm) 




1/3 

(0,0,0) o 





A PROBABILISTIC REGULATORY NETWORK FOR THE HUMAN IMMUNE SYSTEM. 3 

In order to study the dynamic we need the Transition Matrix of the system, 
because we have two functions acting on the set of states. So the dynamic will be 
study using Markov Chains. First we use the following order for the states, but 
this is not the only possibility 

(0,0,0) (0,0,1) (1,0,0) (1,0,1) (0,1,1) (1,1,0) (0,1,0) (1,1,1) 
1 2 3 4 5 6 7 8 

The matrix is constructed in the following way a(i.j) is the probability to have 
the arrow that it is going from i to j, then a( 2j 3) = p{f2) = 1/3, but a( 3j3 ) = 
p{h) +p(fi) = 1 because the two functions are going from (1,0,0) to (1,0,0). 
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The dynamic of the systems is going to stationary states, or the knowing by the 
equilibrium of the system, we have this information with the iteration of the matrix, 
that is computing the power of the transition matrix until we have the same vector 
in each arrow of the matrix. In this very simple case we have two separated spaces, 
so our matrix works with two sub matrices of T: Tu and T22, in fact 

10" 

10 

10 

1 _ 

meanwhile for the other submatrix T 22 we have: 

1 " 
10 
10 
1 _ 

This induce that the equilibrium of the systems is the following 

7T = (0,0, .5, 0,0, .25,0, .25) 

That is, the system is going to the boolean vector (1, 0, 0) with a probability of .5, 
to the boolean vectors (1, 1, 0), (1, 1, 1), with probability of .25. We consider that, 
additionally some part of this systems is going from (1, 1,0) to (1, 1, 1), and from 
(1,1,1) to (1,1,0) continuously. 

3. Projection and Reduction of networks 

Reduction of a network is our interest in this section, we use a bigger net than 
the one in the last section. In particular the mathematical concept that permit us 
to do that is the called homomorphism. In particular here we present a projection, 
that it is an homomorphism which reduce the network to smaller and have the 
very important property to have the similar state of equilibrium, for mathematical 
background see [U [2] . For other approach of the concept of homomorphism, that 





4 



MARIA A. AVINO-DIAZ 




Figure 1. Projection 

had applications to dynamical networks, see [H El ?]. We have the last example 
and now, the following network is defined as follows: A, 

2 • .4 

A \ | 
f • .3 

Each vertex has two values, that is our network is boolean Z 2 = {0, 1}, then the 
functions act on the set Z 2 4 , that has 16 elements or states. The functions are the 
following: 

g 1 (x 1 ,x 2 ,x 3l x 4 ) = (x 1 ,x 2 ,x 1 x 2 ,x 2 xZ), 52(^1, x 2 , x 3 , x 4 ) = (l,x 2 ,xix 2 ,x 2 xi) 

33(0:1, X 2 , X 3 , X 4 ) = (xi,X 2 ,X 2 ,X 2 X^), 94,(xi, X 2 , X3, Xi) = (l,X 2 ,X 2 ,X 2 Xl) 

The homomorphism is the following: 

h(xi,X 2 ,X 3 ,X4) = (xi,x 2 ,x i ) 

This function satisfies the following properties, that are called commutative dia- 
grams: 

£ 2 > £ 2 L 2 > /L 2 

z 2 3 h > z 2 3 z 2 3 h > z 2 3 

f7?4 92 ™ 4 to 4 94 ™ 4 
Z 2 > £ 2 £ 2 > IL 2 



A y h { h { 



z 2 3 12 > z 2 3 z 2 3 12 > z 2 3 

So, h is a structural homomorphism, or an homomorphism of PBN. The proba- 
bilities are p(gi) = p(g 3 ) = 2/3, and p(g 2 ) — p{g 4 ) — 1/6, of course the sum 
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is 1 but this is our suggestion, because we can have others probabilities. In 
fact, using the Theorem in Section 4 [2], this network is going to similar states, 
that is the boolean vectors (0,0,0,0), (0,0,1,0), (1,0,0,0), (1,0,1,0), (0,0,0,1), 
(0, 0, 1, 1), (1, 0, 0, 1), (1, 0, 1, 1) are going to the state (1, 0, 0, 0), and the boolean 
vectors (0,1, 1,0), (1,1,0,0), (0,1,0,0), (1,1,1,0), (0,1,1,1), (1,1,0,1), (0,1,0,1), 
(1, 1, 1, 1) are going to the states (1, 1, 1,0), and (1, 1, 1, 1), and the equilibrium is 
obtained after several iterations of the functions, that is the powers of the transition 
matrix give the following information. 

7t = (0, 0, .5, 0, 0, 0, 0, 0, 0, 0, .25, 0, 0, 0, 0, .25) 

4. The immune system: modeling the B/T-cells interactions 

The following model appears in [9], under the information that it is the model 
B/T — cells, of Kaufman, Urbain, and Thomas, described in 1985. Here we intro- 
duce an important applications of the model PBN to the understanding the immune 
system. 

Antibody = a ^ — > 7), — cells O 
Model B/T - cells \ J, / \\ 

B — cells^ T s — cells O 

The original model had a complex presentation, using a parameter e, for the antigen, 
that has only two values, that is antigen is assumed to be either present or absent 
and hence is represented by e, a binary parameter. The binary variables are b, 
the B-cell population, a, the secreted antibody concentration, and h and s, helper 
and suppressor T-cell populations, respectively. We use for boolean functions, the 
polynomial representations over the field Z2, then we have the following functions: 

fb(b, h, s, a) — eh; fh(b, h,s,a) = es + h + esh; 

f s (b, h, s, a) = h + s + hs; f a (b, h, s, a) = ebh. 

that had obtained by the digraph in the above diagram, called " Model B/T-cells". 
Because the parameter e takes two values 0, and 1, we obtain two functions: 

fo(b, h, s, a) = (0, h,h + s + hs, 0); 

fi (b, h,s,a) = (h,s + h + sh, h + s + hs, bh) 

Then, our space has 2 4 = 16, and we can see how the system is moving to the 
equilibrium of the system It is easy to see that there are three steady states for 
the function / , that is when the antigen is absent, and two steady states for the 
function fi, that is when the antigen is present. In [5], it is considered the function 
/o working and the same time of the function fx but they do not use probabilities, 
so the biological conclusion is not supported by a good description of the activity 
in the net. They consider that the stationary points of /o are the virgin states, that 
is they do not have memory about the antigen, because in those states there are 
only helpers and suppressors cells. But with our analysis we obtain a new results 
Using Markov Chains, and assigning the same probability to each function, that 
is .5, the equilibrium of the system is the following 



TT = (.375, 0, 0, 0, 0, 0, 0, 0, 0, 0, .375, 0, .125, .125, 0, 0), 
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f (b,h,s,a) = (Q,h,h + s + hsJJ) 



f^b^h, s,a) — (h, s +h + sh,h + s + hs,bh) 



Figure 2. Espacio de estados del Modelo Booleano B/T-cell 

where the order in the set of states is the following {(0, 0, 1, 0), (1, 0, 1, 1), (0, 0, 1, 1), 
(1,0, 1,0), (0,0,0,0), (0,0,0,1), (0,1, 0,0), (1,0,0,0), (1,0,0,1), (0,1, 1,1), (0,1, 1,0), 
(0,1, 0,1), (1,1, 1,0), (1,1, 1,1), (1,1, 0,0), (1,1, 0,1)} We suggest the following in- 
terpretation of this equilibrium. The states (0, 0, 1, 0), and (0, 1, 1, 0) are the steady 
states when we do not have antigen, and they have more probability to arrive for the 
network. Meanwhile the others two states are (1, 1, 1,0), and (1, 1, 1, 1), when the 
interaction has a complete action in the network, and they had less probabilities, 
because the life in general, is going to the others states, maybe this happen when 
the helper and the suppressors cells are acting but they are not enough strong to 
destroy the antigen. So a better and more complex dynamical system, can work, if 
we consider three possibilities for the antigen, and for all variables in this particular 
model of interaction. 
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